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Abstract 

The NEMO prototype detector, called "NEMO Phase- 1", has been successfully operated at 2000 m depth from 
December 2006 to May 2007. The apparatus comprises a Junction Box and a Mini-Tower hosting 16 optical sensors. 
Preliminary results are presented. Positions of the optical sensors in the Mini- Tower were reconstructed through 
the acoustic positioning system with high level accuracy. Environmental parameters were analyzed. From data 
corresponding to a livetime of 11.3 hours, atmospheric muon tracks have been reconstructed and their angular 
distributions were measured and compared with Monte Carlo simulations. 
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1. NEMO Phase-1 

The NEMO Phase-1 experiment is a technological 
demonstrator for a km 3 underwater neutrino detec- 
tor. Its final aim is to validate key technologies (me- 
chanics, electronics, data transmission, power dis- 
tribution, acoustic positioning and time calibration 
systems) proposed for the km 3 detector. The appa- 
ratus includes prototypes of the critical elements of 
the km 3 detector: the Junction Box (JB) and a pro- 
totype detection unit called Mini-Tower. On 18 t?l 
December 2006 it was successfully installed 25 km 
East off-shore Catania at 2080 m depth at the un- 
derwater Test-Site of the Laboratori Nazionali del 
Sud in Catania. The Junction Box was connected 
to shore via a 28 km electro-optical cable reaching 
a shore station located inside the port of Catania 
where the power supply, the instrumentation con- 
trol and the data acquisition systems are located. 
The Mini- Tower and the Junction Box were inter- 
connected with an electro-optical jumper equipped 
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with wet mateable hybrid connectors. Connection 
operations were performed with an underwater Re- 
motely Operated Vehicle (ROV). The detector was 
turned off on 18" 1 May 2007. For more details see 
the references ([l|; 0) ■ 

In this work, preliminary results from the NEMO 
Phase-1 detector are presented. Data from the 
acoustic positioning system and environmental pa- 
rameters were analyzed. Atmospheric muon tracks 
have been reconstructed and results were compared 
with Monte Carlo simulations. 

2. The Mini- Tower 

The prototype detection unit built for NEMO 
Phase-1 is a "Mini-Tower" 0) of 4 floors; each floor 
is made of a 15 m long structure hosting two opti- 
cal modules (one down-looking and one horizontal- 
looking) at each end (4 OM per floor). The floors are 
vertically spaced by 40 m. Each floor is connected 
to the following one by means of four ropes that are 
fastened in a way that forces each floor to take an 
orientation perpendicular with respect to the adja- 
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cent (top and bottom) ones. An additional spacing 
of 80 m is added between the tower base and the 
lowermost floor to allow for a sufficient water vol- 
ume below the detector. 

In addition to the 16 Optical Modules (OMs) 
the instrumentation installed on the Mini- Tower 
includes several sensors for calibration and envi- 
ronmental monitoring. Two hydrophones (labeled 
HO and HI) are installed at the end-points of each 
tower-floor, close to the optical modules, and two 
additional monitor hydrophones (called "monitor- 
ing station" ) are installed on the tower base. The hy- 
drophones, together with an acoustic beacon placed 
on the tower base, and other three beacons installed 
on the sea bed, are used for precise determination of 
the tower position by means of time delay measure- 
ments of acoustic signals. The other environmen- 
tal instruments are: a Conductivity- Temperature- 
Depth (CTD) probe used for monitoring the sea- 
water temperature and salinity, a light attenuation 
meter (C-Star) and an Acoustic Dopplcr Current 
Profiler (ADCP) that provides continuous monitor- 
ing of the deep sea currents along the whole tower 
height. In addition, a compass-tiltmeter board is 
placed on each floor in order to measure the tower 
floor inclination and orientation. 

3. Acoustic Positioning System 

A mandatory requirement for the muon tracking 
is the knowledge of the absolute position of the op- 
tical sensors. In fact, while the position and orien- 
tation of the tower base is fixed and known from its 
installation, the rest of the structure can bend un- 
der the influence of the sea currents. A precise de- 
termination of the PMTs position is obtained with 
a system based on time of flight measurements of 
acoustic pulses between a fixed sea-floor array of 
acoustic beacons (that form a Long Base Line, LBL) 
and the target hydrophones located on the detector 
structure. Hydrophone positions are deduced apply- 
ing a triangulation method. Measurements were per- 
formed at regular intervals of ~1 second. Time of 
flight (TOF) is defined as the difference between the 
time of arrival (TOA) of the acoustic signal on the 
hydrophone and the time of emission on the beacon 
(TOE). 

To recognize the beacon pulses a technique called 
Time Spectral Spread Codes (TSSC) was adopted 
in this application. Each beacon transmits a pattern 
of 6 pseudo-random pulses (spaced by ~1 sec) that 
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Fig. 1. Distance H0-H1 measured for floor 2. Each point 
represents the average distance over a period of 300 s. For 
comparison the distance measured during the integration 
phase (d = 14.25 ± 0.01 m, solid line) is shown. 

is different from the others. Each pulse length was 5 
ms and the sequence of pulses was built in such a way 
to avoid overlap between two consecutive pulses. In 
this way a typical beacon pulse sequence can be rec- 
ognized without ambiguity and all the beacons can 
transmit their characteristic pulse sequence at the 
same acoustic frequency. This represents a real ad- 
vantage since all beacons can be identical except for 
the software configuration that defines the pulse se- 
quence, and receivers can be sensitive to only one 
acoustic channel. 

In order to obtain the required accuracy of 15 cm 
(comparable with the dimensions of the PMT) the 
time of flight must be evaluated with an accuracy of 
the order of 10~ 4 sec. To achieve this goal an accu- 
rate calibration of the LBL must be performed and 
the clock drift of the stand-alone beacons must be 
taken into account. The beacon absolute positions 
and relative distances were determined, acoustically, 
at the time of detector installation, using a ROV 
equipped with a 32 kHz pinger, GPS time synchro- 
nized, and a high accuracy pressure sensor. Once 
the distances between the beacons and the monitor- 
ing station, and the sound velocity have been de- 
termined, one is able calculate the TOF. The TOE 
of the beacon pulse, in the Master Clock reference 
time, is obtained measuring the TOA of this pulse at 
the monitoring station. This procedure allows also 
to compensate the clock drift of the stand-alone bea- 
cons (about 20 ns/s) during the whole live time of 
the apparatus. Further, in order to be able to merge 
in post-processing positioning data together with 
optical modules detection information, both were 
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time stamped with a universal time reference tag 
(the Universal Time Coordinate, UTC). 

The relative distances between beacons were mea- 
sured with an accuracy of ~ 10 cm. In order to es- 
timate the accuracy of the positioning system, dis- 
tances between hydrophones HO and HI on the floor 
were measured. In Fig[TJ as an example, the distance 
H0-H1 measured for floor 2 is shown. Each point 
represents the average distance measured over a pe- 
riod of 300 s, in the time interval from I s * February 
h.17 to I s * February h.23 (6 hours). The mean value 
of the measured distance is 14.24 ± 0.06 m (dashed 
line). This value is compared with the distance H0- 
Hl, equal to 14.25 ± 0.01 m (solid line), measured 
on-shore, during tower integration. This result indi- 
cates that the obtained accuracy in the determina- 
tion of hydrophones positions is less than 10 cm, a 
value better than the design requirement (15 cm). 

Data of acoustic positioning system, for the full 
period of the detector operation (from 18 th Decem- 
ber 2006 to 18 th May 2007), were extensively ana- 
lyzed. Tower positions were reconstructed and move- 
ments, as a function of time, were studied within a 
long and short time scale (|4|). 

4. Underwater sea currents 

The Acoustic Doppler Current Profiler (ADCP) 
provides continuous monitoring of the deep sea cur- 
rents. This instrument, mounted on the 4 th floor, 
downward oriented, measures the underwater cur- 
rents below that depth, in a range of about 150 m. 
The ADCP data in the period of the detector data 
taking were analyzed. 

In Fig[2] the magnitude and direction of the sea 
current, measured with the ADCP, are shown. In 
this plot are reported the values averaged in 5 bin of 
depths, corresponding to depth values between 12.5 
and 25 m below floor 4. The recorded current mag- 
nitude is always less than 10 cm/sec and the average 
direction is ~ 180°, this means that the current is 
directed towards the South. 

5. PMTs background counting rate 

Each of the 16 OM hosts a 10" PMT Hamamatsu 
R7081-Sel. The instantaneous rate value is calcu- 
lated by the Front-End board of the PMT averag- 
ing, in a time window of 1 /is, all the hits whose am- 
plitude exceeds a given threshold equivalent to 0.3 
(spe). Background hit rates, for the full period of 
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Fig. 2. Magnitude and direction of the underwater sea cur- 
rents, as a function of time, measured by the ADCP. 
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Fig. 3. Histogram of photon hit rate distribution for a PMT 
located on the 4 th floor, recorded on W th -20 th January 2007. 

the detector operation, were analyzed. The average 
rate of photon hits on all PMTs was measured to be 
ss 80 kHz as expected from 40 K decay plus a contri- 
bution due to diffuse bioluminescence. The fraction 
of time in which a PMT records a high photon rate, 
due to bioluminescence bursts, is typically of the or- 
der of few percent, which implies a negligible dead 
time for the detector in the search for muon events. 

In Fig[3l as an example, the histogram of the rate 
distribution is shown for a PMT located on 4 th floor, 
in the time interval between 10 t?l -20" 1 January 2007. 
In order to disentangle the slowly variable contribu- 
tion to optical background due to 40 K decay and dif- 
fuse bioluminescence, the so called baseline, defined 
as the median value of the rate distribution, was 
calculated. The baseline of the PMT rate was 73.5 
kHz. For comparison the mean value of the distri- 
bution and the average value of the Gaussian func- 
tion fitting the peak are reported. This distribution 
shows also a tail extended to several hundreds kHz 
due to intense bioluminescence bursts. To evaluate 
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Fig. 4. Top panel: PMT baseline rate as a function of time, 
evaluated in term of the median, mean of photo hit rate 
distribution and as average value of the Gaussian fit of the 
histogram peak (see Fig[3j) - Bottom panel: burst fraction 
evaluated as the percentage of time in which the rate is larger 
than 1.2 times the median rate value or larger than 200 kHz. 

the contribution of the biohiminescencc bursts, the 
so called burst fraction, defined as the percentage of 
time for which the rate exceeds a given threshold, 
was calculated. Burst fraction was calculated in two 
different ways: as the percentage of time in which the 
rate exceeds 200 kHz and as the percentage of time 
in which the rate exceeds 1.2 times the median rate 
value. For the data shown in FigGH the percentage of 
bursts larger than 200 kHz is 0.3% while the percent- 
age of bursts larger than 1.2 times the median rate 
value is 2.6%. In FigHthe PMT rate baseline and 
the burst fraction as a function of time, calculated 
averaging in a time interval of 30 min are shown. 



6. Atmospheric muon track reconstruction 

Data recorded during 23 rd and 24 th January 2007, 
corresponding to a livetime of 11.3 hours, were ana- 
lyzed. From the data set analyzed, 2260 atmospheric 
muon events were reconstructed, corresponding to 
a reconstruction rate of 0.056 Hz, and their angular 
distribution was measured (0). 

For comparison, a Monte Carlo simulation of the 
response of the detector to atmospheric muons was 
carried out. A total number of 4 x 10 7 atmospheric 
muon events were simulated with MuPage (j5j), cor- 
responding to a livetime of 4.08 x 10 4 s (~ 11.3 
hours). The energy was sampled between 20 GeV 
and 500 TeV, while the zenithal angle was gener- 
ated in the range [0°, 85°]. The real light absorption 
length spectrum measured at the Test-Site, the op- 
tical background evaluated from the measured PMT 



data rate and the NEMO Phase- 1 DAQ electron- 
ics read-out, were taken as input for the simulation. 
The geometry of the detector was simulated averag- 
ing, the PMTs positions reconstructed in the whole 
day in which the data run was acquired, using the 
acoustic positioning system data. The Test-Site on- 
line trigger settings, used during the data acquisi- 
tion, were also simulated. In order to reject the back- 
ground hits in the whole sample of hits recorded 
in an event, off-line trigger and event selection cuts 
were applied, both in the simulated data and real 
data analysis. Muon tracks reconstruction was per- 
formed using the reco code (@) . 
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Fig. 5. Distribution of reconstructed track directions as a 
function of the cosine of the Zenith angle (cos0^ ec ) (top 
panel) and of the Azimuth angle (0^ ec ) (bottom panel). The 
used likelihood quality cut is equal to -10. 



The angular distributions of reconstructed tracks 
are shown in Fig[5] as a function of the cosine of 
the Zenith angle (cos 6^ ec ) and of the Azimuth angle 
(0^ ec ). The data are shown in comparison with sim- 
ulations. The plotted events were selected applying 
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a quality cut in the likelihood distribution (L) for 
log 10 (L) /Nuoj7>-10, both for data and simulations. 

Despite only a small fraction of the data acquired 
was analyzed in this work, the comparison between 
real data and simulations shows excellent agree- 
ment. The analysis of the whole data set acquired 
by the NEMO Phase- 1 detector (corresponding to 
a live time of ~ 3 months) is in progress. 

7. Conclusions 

The NEMO Phase- 1 detector was operated from 
18 th December 2006 to 18 th May 2007. The posi- 
tions of the hydrophones in the tower were calcu- 
lated with an accuracy smaller than 10 cm. The mag- 
nitude of the underwater current is always smaller 
than 10 cm/sec and the average direction is ~ 180°. 
The baseline for each PMT was measured to be f» 80 
kHz as expected from K decay plus a small contri- 
bution of diffuse bioluminescence. For the data set 
recorded during 23 rd and 24 th January 2007, 2260 
atmospheric muon events were reconstructed. The 
angular distributions of the reconstructed events are 
compared with simulations and a good agreement is 
observed. 
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